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ABSTRACT 

Type-la supernova (SN-Ia) explosions are thought to result from a thermonuclear 
runaway in carbon-oxygen white dwarfs (WDs) that approach the Chandrasekhar 
limit, either through accretion from a companion or a merger with another WD. 
However, it is unknown which of these channels operates in reality, and what are the 
details of the process. I compile from the literature observational estimates of the 
fraction 77 of intermediate-mass stars that eventually explode as SNe-Ia, supplement 
them with several new estimates, and compare them self-consistently. The estimates 
are based on five different methods, each utilising some observable related to the SN- 
Ia rate, combined with assumptions regarding the initial mass function (IMF): the 
ratio of SN-Ia to core-collapse explosions in star-forming galaxies; the SN-Ia rate per 
unit star-formation rate; the SN-Ia rate per unit stellar mass; the iron to stellar mass 
ratio in galaxy clusters; and the abundance ratios in galaxy clusters. The five methods 
indicate that a fraction in the range 77 « 2 — 40% of all stars with initial masses of 
3 — 8Mq (the progenitors of the WDs generally thought capable of exploding through 
the above scenarios) explode as SNe-Ia. A fraction of 77 « 15% is consistent with all five 
methods for a range of plausible IMFs. Considering also the binarity fraction among 
such stars, the fraction of binaries with the appropriate mass ratios, the fraction in 
close initial orbits, and duplicity (every binary can produce only one SN-Ia explosion), 
this implies that nearly every intermediate-mass close binary ends up as a SN-Ia, or 
possibly more SNe-Ia than progenitor systems. Theoretically expected fractions are 
generally one to two orders of magnitude lower. The problem could be solved: if all 
the observational estimates are in error; or with a "middle-heavy" IMF; or by some 
mechanism that strongly enhances the efficiency of binary evolution toward SN-Ia 
explosion; or by a non-binary origin for SNe-Ia. 
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1 INTRODUCTION 

Type-la supernovae (SNe-Ia) play a central role in astro- 
physics, not only as distance indicators for cosmography, 
but also as major contributors to cosmic metal production 
and distribution. The occurrence of some SN-Ia explosions in 
early-type galaxies, having no trace of a young stellar popu- 
lation, implies that at least some of these events are derived 
from old stars, likely low-mass stars or white dwarfs (WDs). 
Detailed analysis of the optical spectra and light curves of 
SNe-Ia (e.g., Mazzali et al. 2007), and of some nearby SN 
remnants thought to have been SNe-Ia (e.g., Badenes et al. 
2006), are consistent with a scenario in which a C-l-O WD 
approaches or passes the Chandrasekhar mass limit, Afch, 
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initiating carbon ignition under degenerate conditions. This 
leads to thermonuclear runaway and incineration of much of 
the WD into iron-peak elements, among which the radioac- 
tive ones power the optical light curve (see, e.g., Hillebrandt 
& Niemeyer 2000, for a review). 

However, the nature of the process behind the growth 
toward the Chandrasekhar mass is not known. The two lead- 
ing competing scenarios are accretion from a close binary 
companion - the single-degenerate (SD) scenario (Whelan 
& Iben 1973; Nomoto 1982) or merger with another WD, 
following loss of orbital energy and angular momentum 
by emission of gravitational waves - the double-degenerate 
(DD) scenario (Iben & Tutukov 1984; Webbink 1984). Both 
scenarios envision passage of the binary through at least one 
common-envelope phase. The complex physics of this stage 
are a major obstacle in calculations of the binary evolution 
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of SN-Ia progenitors. In the SD case, it is also unknown what 
is the nature of the companion - a main sequence, giant, or 
sub-giant star, and what is the form of the accretion flow 
- Roche lobe overflow or a wind. Once conditions for run- 
away have been reached (and these could be sub-Chandra, 
Chandra, or super-Chandra), there are various possibilities 
for where in the WD ignition occurs, and how the burning 
front propagates - various combinations of subsonic ( "defla- 
gration") and supersonic ("detonation") propagation have 
been considered. A number of authors have computed, for 
the various possible progenitor scenarios, the "delay time 
distributions" for the time between the formation of a stel- 
lar population and the explosion of some of its members 
as SNe-Ia (e.g., Greggio & Renzini 1983; Jorgensen et al. 
1997; Sadat et al. 1998; Madau et al. 1998; Ruiz-Lapuente 
& Canal 1998; Yungelson & Livio 2000; Matteucci & Recchi 
2001; Han & Podsiadlowski 2004; Greggio 2005; Belczyn- 
ski et al. 2005). The actual co-existence of several different 
channels that lead to SNe-Ia is possible. In fact, evidence 
has surfaced recently for two different types of SNe-Ia, in 
young and in old stellar populations (see t|2.1l below). 

A better, observationally guided, understanding of what 
objects explode in SNe-Ia, and how it happens, is clearly 
desirable for the confident use of these objects as standard 
candles, and for a coherent picture of metal enrichment. In 
this paper, I revisit the exploding fraction, a simple param- 
eter that can be derived from observations, and which can 
place constraints on the SN-Ia progenitor population. 



2 FIVE ESTIMATES OF THE EXPLODING 
FRACTION 

Any measurement implying a total SN-Ia number from some 
stellar population, relative to some parameter that permits 
normalizing the initial mass function (IMF) of that popu- 
lation, combined with an assumption about the functional 
form of the IMF, can give a direct estimate of the fraction rj 
of stars with initial masses in the range Mmin — A/max that 
eventually explode as SNe-Ia. I compile here five types of 
such measurements, and the explosion fractions that they 
imply, which are also summarised in Table 1. 

In the DD scenario, a minimum zero-age main-sequence 
(ZAMS) mass above Mmin ~ 3Mq is probably required for 
each of the members of a SN-Ia progenitor binary system 
if, after the mass loss involved in post- main- sequence and 
binary evolution, a total WD mass of at least Mch is to 
be attained (e.g.. Branch et al. 1995; Weidemann 2000; Tu- 
tukov & Yungelson 2002) . The maximum initial mass of the 
primary in the binary, Mmax, equals the minimum initial 
mass for eventual core-collapse, generally thought to be in 
the range (8 — 10)Mq (Timmes et al. 1996). It is possible 
that binaries with one or both masses somewhat outside the 
(3 — 1O)M0 range can also produce SN-Ia explosions. For ex- 
ample, in the model calculations by Tutukov & Yungelson 
(2002), the WDs descended from binaries with primaries of 
ZAMS mass > 5Mq and secondaries of ZAMS mass of only 
2.5Mq can still sometimes merge and surpass the Chan- 
drasekhar limit. Tutukov & Yungelson (2002) also assume a 
relatively high minimum mass for core collapse, II.8M0 , and 
therefore quite massive primaries also contribute to their 
calculated SN-Ia rate. The observationally identified pro- 



Method 


»?(%) 


Reference 








la/CC 


8-15 


Mannucci et al. (2005) 


B 


5 — 7 


Dahlcn et al. (2004) 




8-10 


Harris & Tonry (2006) 




6-18 


Scannapieco & Bildsten (2005) 




3-9 


Scannapieco & Bildsten (2005) 




1 - 1.5 


Sullivan et al. (2006) 




3.8 - 4.3 


Mannucci et al. (2006) 


A 


0.8 - 1.7 


Mannucci et al. (2005) 




1 - 1.5 


Sullivan et al. (2006) 




2-6 


Sharon et al. (2007) 




2 - 3.5 


Mannucci et al. (2007) 


Fc 


11 - 20 


Lin et al. (2003) 


Abund 


14 - 40 


De Plaa et al. (2007) 



Table 1. Observational estimates of rj, the fraction of stars with 
ZAMS mass m = 3 — 8 that explode as SNe-Ia. Column header 
explanations: (1) - la/CC - ratio of the type-la and core-collapse 
SN rates in star-forming galaxies; B - the SN-Ia rate per unit 
star-formation rate; A - the SN-Ia rate per unit mass; Fe — the 
ratio of iron mass to stellar mass in galaxy clusters; Abund - 
the number ratio of type-la and core-collapse SNe that have ex- 
ploded in clusters, based on abundance patterns; (2) - Estimate 
of T), in percent, after converting to the adopted IMF parame- 
ters (see text). (3) - Reference for the original estimate, or the 
observational data used in the prsent estimate (see text). 

genitors of several core-collapse SNe have best-fit masses of 
(8 — 10) Mq, but the model uncertainties permit masses as 
high as (12 - 13)M0 (Maund & Smartt 2005; Li et al. 2006; 
Hendry et al. 2006). Nevertheless, among all the DD mergers 
with total mass > Mch in the Tutukov & Yungelson (2002) 
simulations, 85% originate from binaries in which both mem- 
bers had ZAMS masses between (3 — 8)Mq (L. Yungelson, 
private communication), and therefore this is a useful range 
to consider for the progenitors of the large majority of SNe- 
Ia. 

I will henceforth use the following notations. Stellar 
masses will be expressed in dimensionless form, as m = 
M/Mq. The IMF, dN/dm, is given in terms of stars per unit 
mass interval. I denote with miow and mhi the low-mass and 
high-mass cutoffs, respectively, of the IMF. I will begin by 
assuming the following "standard" parameters: mmin = 3; 
mniax = 8; a Salpeter (1955) IMF, with {dN/dm) oc m"^'^^; 
?Tiiow ~ 0.1; and mhi ~ 100. In i )2.6l I will examine the con- 
sequences of varying these parameters. The treatment will 
first focus on the DD scenario, but in 33] I will investigate 
how the results are affected in the SD case. 

2.1 The local ratio of type- la to core-collapse SN 
rates 

Mannucci et al. (2005) have reanalysed the SN survey data 
described in Capellaro et al. (1999), to obtain core-collapse 
and type-la SN rates for distinct galaxy types. They, and 
several subsequent studies (Scannapieco & Bildsten 2005; 
Mannucci et al. 2006; Sullivan et al. 2006) , have found that 
the SN-Ia rate in star-forming galaxies is dominated by a 
"prompt" progenitor component. This component has a typ- 
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ical delay time between stellar formation and SN-Ia explo- 
sion of ~ 10* yr or less. This leads to a linear dependence 
of SN-Ia rate on the star formation rate of the host galaxy 
population. In star-forming galaxies, Mannucci et al. (2005) 
find a ratio of Ni^/Ncc ~ 0.35 ± 0.08 between the type-la 
and core-collapse rates. This ratio permits a first estimate of 
the fraction of the potential type-la progenitor stellar pop- 
ulation that actually explodes through the prompt type-la 
channel: 



^la/CC 



A^ia ff{dN/dm)dm 



(1) 



Ncc J^{dN/dm)dm 

For the standard parameters, the ratio of the two integrals 
equals 0.332, giving rj = 0.12 ± 0.03. In the integral in the 
numerator, the upper mass limit is 50, since more massive 
stars may collapse directly to black holes without a SN ex- 
plosion. Lowering the upper limit to 25 reduces rj by 14%, 
and raising it to 100 increases it by 9%. Including this un- 
certainty gives, for the standard parameters, an estimate of 
Vi./cc = 0.08 - 0.15. 

2.2 The SN-Ia rate per unit star-formation rate 

A "prompt" SN rate per unit star formation rate, B (adopt- 
ing the notation introduced by Scannapieco & Bildsten 
2005), has dimensions of a SN number divided by a stel- 
lar mass. The value of B equals the total number of prompt 
SNe-Ia that have exploded among an actively star-forming 
population since star formation began, divided by the total 
stellar mass that has formed. Assuming that the population 
is being observed at least > 10* yr after the onset of star for- 
mation, the fraction of stars with initial masses in the range 
w-min — w-max that eventually explode as SNe-Ia is therefore 



riB = B 



J^^^^ m{dN / dm)dm 
f^{dN/dm)dm 



(2) 



For the standard parameters, the ratio of the two integrals 
equals 47.6. 

It is important to remember that SN rate measurements 
are not normalized by a directly measured star formation 
rate. Rather, the luminosity and the color, or set of col- 
ors, of the host galaxy population are used to derive a star 
formation rate, usually by comparison to spectral synthe- 
sis models. The luminosity and colors of a stellar popula- 
tion at a given time is dominated by the most massive stars 
that are evolving off the main sequence at that time, while 
the mass of the population is dominated by the least mas- 
sive, but most numerous, stars. The derivation of a star for- 
mation rate therefore requires the assumption of an IMF 
for the largely unobserved low-mass population. For self- 
consistency, the right side of Eq. [2]should then be multiplied 
by a correction factor 



Jr 



m{dN/dm)oTsdm 



yimf ,org 



J^'° m{dN/dm)dm 



(3) 



where mto is the turnoff mass, and the "org" subscript de- 
notes the original IMF assumed in the published rate mea- 
surement. Previous studies that have derived SN rates have 
generally assumed the same, or similar, IMF parameters as 
the standard ones I have chosen, making (?imf,org ~ 1. I will 
note, below, when this is not the case. 



Dahlen et al. (2004) measured the SN-Ia rate vs. red- 
shift z, out to z ^ 1.6. They took a parametric fit to the 
star- formation rate vs. cosmic time, SFR{t), found by Gi- 
avalisco et al. (2004), convolved it with various parametrised 
trial SN-Ia delay time distributions, and fitted the observed 
SN-Ia rate vs. z. The normalization in this fit, (1.0 — 1.3) x 
10"^Mq^ for the various delay functions, gave the number 
of SNe-Ia per unit formed stellar mass. As in Eq. [2] above, 
Dahlen et al. (2004) multiplied this by the ratio of the same 
two integrals (they used mhi — 125, rather than mhi — 100, 
which makes a negligible difference), to obtain rjB = 5% to 
7%. 

Barris & Tonry (2006) likewise compared the SFR{t) 
of Giavalisco et al. (2004) to the SN-Ia rates they measured 
themselves, but used an interpolated form of the SFR{t), 
rather than a fitted functional form, as done by Dahlen et 
al. (2004). Furthermore, they limited their comparison to 
z < 1, and rather than convolving the SFR{t) with a de- 
lay function, they shifted it by a set delay (equivalent to 
convolution with a shifted 5-function, which is reasonable 
for the short delays they considered). They found a best-fit 
delay of — 1 Gyr, a number of SNe-Ia per stellar mass of 
(1.7 - 2.2) X 10-^M-\ and rjB = 8% to 10%. 

Scanappieco and Bildsten (2005) used the Dahlen et al. 

(2004) measurements at z ^ 1 of type-la and core-collapse 
SN rates, and the Giavalisco et al. (2004) SFR{t), to derive 
B = (2.6 ± 1.1) X lO~^SNeM0"\ For the standard IMF 
parameters, this gives rjB = 12 i: 6%. Alternatively, Scanap- 
pieco and Bildsten (2005) multiplied the SN-Ia rate mea- 
sured by Mannucci et al. (2005) in starburst galaxies by the 
age of the stellar population of these galaxies (as estimated 
based on their colors), to obtain B = (1.2±0.6) x 1O"^M0"\ 
which I translate to rys — 3% to 9%. 

Sullivan et al. (2006) estimated star-formation rates, 
based on galaxy broad-band spectral energy distributions, 
for the SN host galaxies in the Supernova Legacy Survey, 
at 0.2 < z < 0.75. They found a prompt component of the 
rate of B = (3.9 ± 0.7) x 10"^ Mq''^. To obtain their mass 
estimates, Sullivan et al. (2006) assumed a Kroupa (2001) 
IMF, which has a slope of a = —2.3 above a mass of m = 0.5 
and a — —1.3 below it. If the Kroupa (2001) and Salpeter 
(1955) IMFs are both normalized at a mass of 8 (which cor- 
responds roughly to the B-type stars that dominate the light 
of young stellar populations) , the ratios of the respective to- 
tal masses are 0.69. Converting from the Kroupa (2001) IMF 
to the standard IMF parameters therefore increases the to- 
tal stellar mass formed, and hence reduces B by 0.69. The 
measurement of B by Sullivan et al. (2006), adapted to the 
standard IMF parameters, then implies = 1 — 1.5% 

Mannucci et al. (2006) have tested a large variety of de- 
lay time distributions, some parametrized and some phys- 
ically motivated, to fit simultaneously: (a) the SN-Ia rate 
vs. redshift (by convolving, again, the parametric fit to the 
SFR{t), found by Giavalisco et al. 2004, with the delay 
function); (b) the SN-Ia rate per unit mass as a function 
of galaxy colors measured by Mannucci et al. (2005); and 
(c) the dependence of SN-Ia rate on galaxy radio loudness 
found by Delia Valle & Panagia (2003) and Delia Valle et al. 

(2005) . As before, the normalization of the fit to the SN-Ia 
rate vs. redshift gives rj. For both of their best-fitting mod- 
els - a physical SD model by Belczynnski et al. (2005), and 
a parametrized two-component delay function, consisting of 
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equal contributions from a prompt component and from a 
delayed component having an exponential decay with char- 
acteristic time 3 Gyr - they found values of 77 « 4%. 

Thus, for the standard parameters, these six estimates 
imply TjB in the range 1 — 10%. This can be compared to 
Vin/cc = 8 — 15%, found above, as both relate to the same 
prompt SN-Ia component. 



2.3 The SN-Ia rate per unit stellar mass 

A third estimate of rj can be obtained from the "tardy" com- 
ponent of the SN-Ia rate identified by Mannucci et al. (2005), 
and which dominates in old stellar populations, whose ages 
are ~ 10 Gyr. This is effectively a "DC" component in the 
SN rate, whose level is proportional to the total stellar mass, 
and is denoted by A (again following Scannapieco & Bild- 
sten 2005). Assuming that the maximum delay time between 
star formation and SN-Ia explosion for this component is 
~ 10 Gyr, then A x 10 Gyr is the number of SNe-Ia per 
unit stellar mass that explode through this channel, and the 
SN-Ia explosion fraction in this case is 



in L , midN I drn)dm 



/„ {dN / dm)dm 



(4) 



Formally, the integral in the numerator should run up 
only to the main-sequence turnoff mass in old stellar popu- 
lations, m ~ 1. Above m ~ 1, the numerator should include 
only the contributions from stellar remnants - WDs, neu- 
tron stars, and black holes. More to the point, like Eq. [2j 
Eq. |4] needs to be corrected by the factor pimf.org (Eq. [3]) 
appropriate for the IMF assumed in the original SN-Ia rate- 
per-unit-mass calculation, when deriving the mass of the 
stellar population from its observed luminosity and colors. 
Often, the total formed mass, rather than the mass in stars 
and in remnants, has been used. 

Note that Eq. [4] assumes that, ~ 10 Gyrs after the for- 
mation of the stellar population, we are seeing the last SNe- 
Ia that form through the tardy channel. If the SNe, in reality, 
continue to explode at a constant rate for, say, 20 Gyr (i.e., 
the supply of SN-Ia progenitors has been only half used up 
after 10 Gyr), then ri would be doubled. Alternatively, if the 
tardy component is not constant, but, say, decays exponen- 
tially with a timescale of 10 Gyrs, then Eq. |3]is exact. 

Mannucci et al. (2005) derived a local value of the 
mass-normalized SN-Ia rate in early-type galaxies of ^ = 
(3.8 ± 1.3) X 10"''yr-i(10^° Mq)-^ (after converting their 
measurement from Hq of 75 to 70 km s~^Mpc~^). Their 
mass determination was based on the calculations of Bell 
& de Jong (2001), who defined relations between the mass- 
to-light ratios and the optical-to-near-IR colors of synthetic 
galaxies. Bell & de Jong (2001) assumed a Salpeter (1955) 
IMF, but scaled down by a factor 0.7, which they denoted 
a "diet" Salpeter IMF. Including this factor in Eq. 3] we 
obtain tja = 0.8 - 1.7% 

Sullivan et al. (2006) found, based on the early-type 
galaxies in the Supernova Legacy Survey, a tardy component 
of the SN-Ia rate of ^ = (5.3 ±1.1) x lO'^yr"^ (lO^^Mo)-^ 
As noted above, in obtaining their mass estimates from stel- 
lar synthesis modeling, they assumed a Kroupa (2001) IMF, 
and only the mass in surviving stars. If the Salpeter (1955) 
and Kroupa (2001) IMFs are both normalized at a mass 



of 1 (which corresponds roughly to stars that dominate the 
light of old stellar populations), the ratio of the total stellar 
masses below m = 1 equals 2.0. Converting from the Kroupa 
(2001) IMF to the standard IMF parameters therefore in- 
creases the total stellar mass formed (and hence reduces A) 
by a factor of 2, giving 77A = 1 — 1.5%. 

Sharon et al. (2007) have recently measured a SN-Ia rate 
in massive clusters at z « 0.1 (based on the SN survey by 
Gal- Yam et al. 2007), oi A = 9.8t| j x lO-''yr-i(lO^°M0)-\ 
Since there is little star formation in clusters at these red- 
shifts, this is effectively another measurement of the tardy 
A component of the SN-Ia rate. Following Mannucci et al. 
(2005), the mass determination was based on the mass-to- 
light ratios found by Bell et al. (2003), who assumed a 
Salpeter (1955) IMF, scaled down by a factor 0.7. Including 
this factor in Eq. UJ we obtain r/A = 2 — 6% 

Mannucci et al. (2007), again using the Cappellaro et al. 
(1999) sample, find a rate in local early-type cluster galaxies 
of A = e.ela o X lO"''yr"^(lO^°M0)-^ corresponding to 
r/A = 2 - 3.5%' 

We thus obtain values of r]A between 1% and 6% for the 
tardy component. The total fraction of the initial m = 3 — 8 
stars in a stellar population that explode as SNe-Ia through 
either one of the prompt or tardy channels is 7? = 7?a + »7s or 
f) — »7A+T;ia/cc- Taking the extreme values found above, this 
gives a range for rj of 2% to 21% for the standard parameters. 



2.4 The iron to stellar mass ratio in galaxy 
clusters 

The X-ray spectra of massive galaxy clusters reveal a large 
mass of iron, most of it in the hot intra-cluster medium 
(ICM, e.g., Balestra et al. 2007, and references therein). A 
tight correlation is observed between the iron mass and the 
total luminosities of the early-type galaxies in clusters, while 
no correlation is observed with the sub-dominant late-type 
galaxy luminosity. This strongly suggests that the source 
of the iron is SNe from the stellar population whose low- 
mass component currently constitutes the early-type galax- 
ies (rather than, say, the late-type galaxies, or infall of exter- 
nal material into the cluster). Assuming the standard IMF, 
it can be shown (e.g., Maoz & Gal- Yam 2004) that core- 
collapse SNe could have produced only about 20% of the to- 
tal mass of iron in the ICM and in the galaxies themselves. 
If the remaining 80% of the iron was formed by SNe-Ia, then 
the ratio of this iron mass to the current stellar mass gives 
another, independent, estimate of 77: 



0.8Mfc fg°° ■m{dN/dm)dm 
Yf,M, J^{dN/dm)dm ' 



(5) 



where Mpc/Mt, is the ratio of total iron mass to current total 
stellar mass, and Ypo is the mean iron yield of a single SN- 
Ia. The first term in the product on the right side is just 
the cumulative number of SNe-Ia that have exploded in the 
cluster, divided by the stellar mass. 

From analysis of X-ray, optical, and IR measurements, 
Lin, Mohr, & Stanford (2003) have estimated the iron-to- 
stellar mass ratio of clusters. Their figure 9, with their as- 
sumed Solar iron abundance of Zq = 1.814 x 10"^ (An- 
ders & Grevese 1989), suggests Mfc/M, = (3 ± 0.5)^0 = 
(5.4 ±0.9) X 10"^. Lin et al. (2003) found the stellar masses 
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in their sample clusters based on the 2MASS survey Ji'-band 
luminosities of the cluster galaxies. The spiral galaxy stellar 
mass-to-light ratio, M/L, was obtained from the relations 
with luminosity and color of Bell & de Jong (2001), who 
assumed the "diet" Salpeter IMF (i.e., scaled down by a 
factor 0.7, see above). For the ellipticals, Lin et al. (2003) 
used the dynamical estimates by Gerhard et al. (2001) of 
M/L as a function of luminosity in a sample of 21 ellipti- 
cals. Lin et al. (2003) found, for the clusters as a whole, a 
range of if-band M/L between 0.7 and 0.8, in solar units. 
They noted the similarity of this value to the mean value of 
0.73 obtained by Cole et al. (2001) for the 2dF galaxy red- 
shift survey. In deriving this value. Cole et al. had assumed 
a Kennicutt (1983) IMF, which has a slope of a — —1.4 for 
m < 1 and q: = — 2.5atm>l. Converting from a Kenni- 
cutt (1983) to a Salpeter (1955) IMF increases the Jf-band 
M/L of an old stellar population by a factor of 2.0 (Bell et 
al. 2003). Thus, for our standard IMF parameters, the Lin 
et al. (2003) iron-to-stellar mass ratio is reduced by a factor 
of 2, to Mfc/M. = (2.7 ± 0.5) x 10"^ 

The mean iron yield of a SN-Ia is observationally con- 
strained by SN-Ia light curves and spectra, as well as from 
SN la model calculations, to be Ypo ~ 0.7 ± O.IM© (e.g., 
Thielemann, Nomoto, & Yokoi 1986; Nugent et al. 1997; 
Contardo, Leibundgut, & Vacca 2000; Mazzali et al. 2007; 
although the iron mass yield among individual SNe-Ia likely 
has a dispersion of a factor of a few, for the present estimate 
we are actually interested only in the mean, and the error 
on the mean). Accounting for the uncertainties then gives a 
range rjpe = 11 — 20%. 

The core-collapse contribution to the iron production 
may have been greater than the 20% stated above, if the 
first stellar generations formed with a top-heavy IMF that 
favored massive stars that exploded as core-collapse SNe, 
or if many of the ancient core-collapse SNe were more effi- 
cient iron producers than present-day core-collapse SNe (see 
Maoz & Gal- Yam 2004; Loewenstein 2006). Either of these 
non-standard scenarios would lower rjFe - On the other hand, 
Kapferer et al. (2007) find, by performing simulated X-ray 
observations of simulated clusters, that X-ray measurements 
systematically underestimate the mass of metals in the ICM. 
Correcting for this effect would raise rjFt^- 

The low SN-Ia rates observed in clusters at z < 1 (Gal- 
Yam et al. 2002; Sharon et al. 2007; Mannucci et al. 2007) 
suggest that only the tardy SN-Ia channel is presently active, 
while much of the iron was produced at early times by a 
prompt channel. The iron mass in clusters is likely a tracer 
of the cumulative products of both the prompt and tardy 
SN-Ia channels. If so, the range of ryFe, 11-20%, should be 
compared to rjA + rjB, or ryx + ?7ia/cCi above, which have a 
range of 2 - 21%. 

2.5 The abundance-based ratio of core-collapse to 
type-la SNe in clusters 

The large effective areas and high X-ray spectral resolu- 
tions of Chandra and Newton-XMM have permitted accu- 
rate measurements of the ICM abundances of additional el- 
ements, beyond iron. De Plaa et al. (2007) have recently de- 
rived abundances for a sample of clusters, and fitted to the 
abundance patterns the theoretical yields of core-collapse 
and type-la SNe. They then derived the relative numbers 



of the two types of SNe, whose cumulative outputs have 
contributed to the observed ICM abundances. Their best 
fit uses, for the SN-Ia contribution, a delayed detonation 
model by Badenes et al. (2003) that Badenes et al. (2006) 
had found provides the best fit to the X-ray emission from 
the Tycho supernova remnant. The best fit to the ICM 
abundances then indicates that the number ratio of the two 
types of SNe that have exploded in clusters is in the range 
Nia/Ncc = 0.5 - 1.20 This leads to a fifth estimate of 77: 

_ Njg C{dN/dm)dm 

which is the same kind of calculation as in Eq. [T] above. 
De Plaa et al. (2007) performed this calculation, based on 
their derived Nja/Ncc ratio, and the same assumed IMF 
parameters, but taking mmin = 1.5 or 0.9 rather than 3, 
JTimax = 8 or 10, and either a Salpeter or a Kroupa (2001) 
IMF. Using mmin = 3, and again allowing also for the un- 
certainty in the upper mass limit for core-collapse SNe, I 
find r?abund ~ 14 — 40%. Like jypc above, ryabund measures the 
exploding fraction through both SN-Ia channels. 

2.6 Dependence of rj on assumed parameters 

For the standard parameters with a Salpeter IMF (i.e., a 
single power law of index a = —2.35), the integrals above 
are trivially solvable, and give the stated fractions. I now 
investigate briefly the effect or varying those parameters, 
within plausible values. 

Most of the forms that have been proposed for the IMF 
differ mainly at m < 1, while at m > 1 they are usu- 
ally single power laws of similar slopes, between —2.2 and 
-2.7 (e.g., Baldry & Glazebrook 2003, a = -2.2; Kroupa 
2001, a = -2.3; Chabrier 2003, a = -2.3; Salpeter 1955, 
a = -2.35; Gould et al. 1997 a = -2.35; Kennicutt 1983, 
a = —2.5; Kroupa et al. 1993, a = —2.7; an exception is 
the Scalo 1998 IMF, which has a = -2.7 for 1 < m < 10 
and a = —2.3 at m > 10). The actual observables in SN 
rate measurements are the ratios of SN numbers of different 
types (la to core-collapse) , or of a SN number and a stellar 
luminosity, the latter dominated either by m « 1 stars in old 
populations, or by m ~ 8 stars in young populations. Nor- 
malizations of SN rates relative to total stellar mass or total 
star formation rate are always obtained by assuming some 
IMF over the entire stellar mass range. However, estimates 
of rj relate only to the stellar populations in the intermediate 
mass (SN-Ia progenitors) and high mass (core-collapse pro- 
genitors) ranges. In my estimates of rj, above, I have taken 
account of the IMFs that were assumed in the original rate 
calculations and converted them consistently to estimates 
based on the standard IMF parameters. Any actual depen- 
dence of 77 on the IMF will therefore be only through the 
form of the IMF at m > 1. 

^ This result is not in conflict with the previous argument, that 
80% of the iron in clusters may come from SNe-Ia. Each SN-Ia 
produces an order-of-magnitude more iron than a core-collapse 
SNe, so both results suggest that SNe-Ia dominate the iron pro- 
duction. However, one must keep in mind that the elemental yield 
of core-collapse SNe based upon explosion models is extremely 
uncertain. 
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Figure 1. Allowed regions in the parameter space of r], the frac- 
tion of all stars with ZAMS mass m = 3 — 8 that explode as 
SNe-Ia, and a, the power-law slope of the IMF at m > 1. The 
two regions shown are based on two methods for estimating the 
fraction of stars that explode only through a "prompt" SN-Ia 
channel, that occurs within ~ 10* yr of star formation, r/g (dotted 
curve) is based on measurement of the SN-Ia rate per unit star- 
formation rate in star-forming galaxies. r/j^/cQ is based on the 
ratio of Type-la and core-collapse SN rates in local star-forming 
galaxies. The vertical extent for each method is based on the 
union of the Icr uncertainty ranges of the different measurements, 
described in the text. The left and right borders encompass the 
range of high-mass slopes of most standard IMFs. 



Figure 2. Same as Fig. 1, but for estimates of rj that reflect the 
total fraction of intermediate-mass stars that explode as SNe-Ia, 
through either a prompt or a tardy channel. The dashed curve, 
labeled Fe, is based on the iron-to-stellar mass ratio in galaxy 
clusters. The dotted curve, for r/abundi is based on the deduced 
ratio of la to core-collapse SN numbers that have contributed to 
observed ICM abundances. The lower border of the solid curve is 
the sum of rig (shown in Fig. 1), and rjj^, the fraction exploding 
through the tardy component, based on SN-Ia rates in old popu- 
lations. The upper border of the solid curve is the sum of »7ia/cCi 
shown in Fig. 1, and tja- A value of ?? fs 15% is consistent with 
all of the methods and most of the IMF slopes. 



IMFs that arc flatter in the mmin— ?Timax range will yield 
either fewer or more progenitors, depending on whether the 
IMF is normalized at the high end or the low end, respec- 
tively. For r)A and r?Fe, the IMF normalization needs to be 
made at m « 1, the stellar mass that dominates the hght 
of the old populations whose mass normalizes those SN-Ia 
number estimates. Therefore, a flatter IMF slope will in- 
crease the number of SN-Ia progenitors and lower rj, and 
vice versa. Conversely, rjB is based on SN-Ia rates normal- 
ized by the light of young stellar populations, dominated by 
the O- and B-type stars at the high end of the SN-Ia progen- 
itor range. A flatter IMF slope will raise ris, and vice versa. 
The estimates of 77ia/cc and r^abund are based on number rar- 
tios of high mass and intermediate-mass stars. A flatter IMF 
slope will always increase the number of high-mass stars rel- 
ative to low-mass stars, and therefore r? will increase in these 
cases. 

Figure 1 shows the ranges of the two estimates of rj for 
the prompt channel, namely ??ia/cc a-nd rjB, and their depen- 
dences on the IMF slope, a, above m = 1. For each method, 
at each value of the IMF slope, the range of permitted rj is 
the union of the la uncertainty ranges of the difi'erent mea- 
surements based on the method. There is a broad range of 
consistency between those two estimates in the range -q — 5 
to 10%, for the various plausible high-mass IMF slopes. Note 
that r/ia/cC) which is proportional to the observed ratio of 
type-la and core-collapse rates, will be systematically bi- 
Eised high if an unaccounted fraction of core-collapse SNe in 
star-forming galaxies are undetected because they explode 
in highly obscured regions. Indeed, a number of studies (e.g., 
MaioUno et al. 2002; Mannucci et al. 2003; Cresci et al. 2007) 
have concluded that a large fraction of core-collapse SNe in 



starburst galaxies are missed by optical surveys. Mannucci 
et al. (2005) show, based on a simple model, that of or- 
der one-half of the local core-collapse SNe could possibly be 
missing. Correcting for this would lower ?;ia/cc by one-half 
and result in a larger overlap in Fig. 1 between the allowed 
areas corresponding to rjE and J7ia/cc- A further reduction 
in 7)ia/cc will result if rn = 8—10 stars form ONcMg WDs, 
rather than core-collapse SNe, which is still the subject of 
theoretical uncertainty (see, e.g., review in Herwig 2005), or 
if low-mass or low-metallicity core-collapse progenitors pro- 
duce faint SNe that are missed by surveys. For example, 
limiting the detectable core-collapse SNe to those with ini- 
tial m > 12 would lower r?ia/cc by another factor of ^ 0.5. 

Figure 2 is the same as Fig. 1, but for the estimates of 
r] that trace the cumulative effect of both the prompt and 
the delayed channels, i.e., tja + Vb, Va + Via/cc, Vve, and 
'Jabund. For claiity, I have plotted the union of the regions 
in parameter space allowed for rjA + Vb (lower limit of the 
region) and rjA + r/ia/cc (upper limit of the region). Interest- 
ingly, a value of ~ 15% is consistent with all of those different 
methods of estimating the total rj, as long as a assumes one 
of its coiivciitioiial values. 

In terms of vaxying the other standard parameters, en- 
larging the range TTT-min ^max will naturally increase the 
number of potential SN-Ia progenitors, and therefore will 
lower ri- Estimates of ri that involve also the core-collapse 
SN rate will be further lowered by the smaller range of core- 
collapse progenitors resulting when mmax is increased. For a 
Salpeter IMF, lowering mmin from 3 to 2.5, 2, or 1.5 lowers r] 
by factors of 0.72, 0.5, or 0.32, respectively. For estimates of 
77 not involving the core-collapse rate, raising mmax to 9, 10, 
or 11 lowers rj by another factor of 0.94, 0.91, or 0.89, respec- 
tively. When the core-collapse rate is involved (?7ia/cc and 
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^yabund), the respGctivG reductions are 0.80, 0.65, and 0.55. 
On the other hand, there is ongoing uncertainty regarding 
the highest initial stellar mass that leads to a CO WD. Re- 
viewing the stellar-evolution models and the observational 
evidence (the initial-final stellar mass relation), Woidcmann 
(2000) argues for a highest initial mass of m = 6 — 7, al- 
though some models (Girardi et al. 2000) find even lower 
limits, of m ~ 5. A shrinking of the m = 3 — 8 progenitor 
range tom = 3 — 5orm = 3 — 6 will naturally raise r/. All 
these changes can thus amount to reductions or increases by 
up to a factor of a few in all of the estimates of rj, derived 
above. 



3 BINARITY, MASS RATIO, SEPARATION, 
AND DUPLICITY 

The fraction r?, discussed above, is the fraction of all stars 
with ZAMS masses in the prescribed range that explode as 
SNe-Ia. To estimate the fraction of all close binary systems 
with members in this mass range that lead to a SN-Ia ex- 
plosion, one most consider the binarity fraction, the binary 
mass-ratio distribution, the separation distribution, and the 
duplicity of the progenitors. 

The binarity fraction among stars is a function of meiss, 
with most low-mass stars being in single systems, while most 
massive stars are members of binaries (see Lada 2006, and 
references therein). Amongst the intermediate-mass popu- 
lation of interest here (B-type and A- type stars), the bi- 
narity fraction (defined as the fraction of systems that are 
binary, rather than single-star, systems) is likely in the range 
b « 1/2 to 1 (e.g., Shatsky & Tokovinin 2002; Kouwenhoven 
et al. 2005; Raines et al. 2006). If a fraction b of systems are 
binaries, and 1 — b are single, the fraction of all stars that is 
in binaries is fa = 2b/{b + 1). Thus, fa « 2/3 to 1. 

Several recent studies have measured the binary mass 
ratio distribution, f{q), where < g < 1 is the ratio of sec- 
ondary to primary mass, for intermediate-mass stars. For 
A-typc stars, Soedorhjolm (2001) finds that f{q) ~ const., 
for separations of 60 — 120 AU. For B-typc stars, Shatsky& 
Tokovinin (2002) find f(q) ~ q"'"'^ to q^""^ at separations 
of 45 - 900 AU. Kouwenhoven et al. (2005) measure q~°-^^ 
at 30 — 1600 AU for A-type and late B-type stars in the Sco 
OB2 association. A distribution f{q) that is flat or slowly 
rising toward small q appears to be generic for a wide range 
of masses and separations. For example, Mazeh et al. (2003) 
find such an f{q) for lower-mass (late G and early K) close 
binaries with periods < 100 days (i.e., separations < 1 AU). 
The fraction of primaries with m = 3 — 8 that have a sec- 
ondary also in this mass range (as required in the DD sce- 
nario) is 

^ l!(dN/dm)dmf^^^,^f{q)dq 
f^idN/dm)dm 

Assuming f{q) oc q^"''' (and normalizing f{q) so that its 
integral from to 1 is unity), we find ft — 0.17. Taking 
f{q) = constant gives ft = 0.29. Thus « 1/6 to 1/3. The 
fraction fi, depends only weakly on the IMF slope in the 
m = 3 — 8 range. 

The initial distribution in physical separation in bi- 
naries, a, is still debated. Oepik (1924) already suggested 



that this distribution has the form /(a) ~ i^^, but other 
forms have been suggested (e.g., Duquennoy & Mayor 1991). 
Poveda et al. (2006) have recently analyzed several sam- 
ples of binaries, and concluded that the Oepik form pro- 
vides an excellent fit over a very wide range, from a = 60 
to 45,000 AU (the upper limit determined for binaries in 
the young Orion Nebula cluster, which have not had time 
to be disrupted by encounters with massive perturbers), 
and suggesting that the same physical process operates at 
all scales. Both of the currently popular scenarios for SNo- 
la, SD and DD, require passage of the progenitor binary 
through at least one common-envelope phase. For initial 
m, = 3 — 8, the maximum radii that single stars achieve 
during their post-main-scqucncc evolutions are 2-5 AU (see, 
e.g. , fig. 1 of Yungolson 2005) . This dictates an initial separa^ 
tion < 2— 5 AU, if the binary is to undergo common-envelope 
evolution. Taking the full range of possible initial separations 
from 2 X lO"'^ AU (contact between two m — 3 ZAMS stars) 
to lO'* '^ AU, and assuming the 1/a distribution, a fraction 
of fc « 0.3 — 0.4 of binaries have an initial separation within 
2-5 AU. If we assume, "generously", that even binaries with 
initial a ~ 50 AU can migrate to smaller separations (e.g., 
via interactions with protoplanetary disks), then fc ~ 1/2. 
Thus, fc is likely in the range 1/4 to 1/2, but probably at 
its low end. 

Finally, every DD binary can lead, at most, to one SN- 
Ia, and therefore the duplicity factor is fd ~ 1/2. 

If rj is the fraction of all intormodiato-mass stars that 
explode as SNe-Ia, the fraction of all intermediate-mass close 
binaries that explode is fia = v/ifafbfcfd)- The range 
formed by the products of the extreme values of fa through 
fd, above, gives fia ~ (12 — 70)77. For all of the observa- 
tionally derived values of rj derived above, a fraction of or- 
der 100% of all intermediate-mass close binaries explode as 
SNe-Ia. For some estimates of r] and of the binary param- 
eters (for example, for 77 — 15%, a value that is consistent 
with all methods, and a product of the binary parameters 
of 70), there are up to an order of magnitude more SNe-Ia 
than potential progenitor systems . 



4 DISCUSSION 

The empirical derivation of a high SN-Ia formation efficiency 
is not new. For example, based on observations by Tam- 
mann (1978) that the type-la and core-collapse SN rates in 
spiral galaxies are equal, Greggio & Renzini (1983) already 
reached the same conclusion. The authors of several of the 
observational studies whose results I have used above have 
themselves made the same point. The novelty in the present 
work is that: (a) I have based my estimate of on a num- 
ber of relatively accurate measurements of extragalactic SN 
rates, and of related observables, that are recently available; 
(b) I have compared among the values of rj and its uncer- 
tainty resulting from difi'erent measurements of the same 
type (e.g., SN la-rate per unit star-formation rate), and to 
those derived by independent methods, all after correcting 
for the different IMF assumptions that went into the vari- 
ous estimates; (c) I have tested for the dependence of the 
results on IMF assumptions, and argued that they depend 
only on the high-mass (m > 1) slope, which is fairly well 
constrained, (d) For all of these estimates, I have incorpo- 
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rated recent results on binarity fraction, binary separation, 
and binary mass-ratio distribution, to estimate the fraction 
of close, intermediate-mass, binaries that explode. 

I have shown that, if our basic ideas about SN-Ia pro- 
genitors are correct, then, at best, the processes leading 
to a SN-Ia explosion are extremely efficient, in the sense 
that most or all of the potential progenitor systems - close, 
intermediate-mass, binaries - do explode as SNe-Ia. At 
worst, some combinations of parameters imply the existence 
of many more explosions than progenitor systems, indicating 
an error in one or more of the assumptions. 



4.1 Theoretical expectations 

What are the theoretical expectations for /ja? Calculations 
in which binaries with a range of initial parameters are fol- 
lowed throughout their stellar and binary evolution, have 
been carried out by, e.g., Iben & Tutukov (1984), Hachisu et 
al. (1996, 1999a,b), Yungelson & Livio (1998, 2000), Langer 
et al. (2000), Tutukov & Yungelson (2002), Hurley et al. 
(2002), Han & Podsiadlowski (2004), Fedorova et al. (2004), 
and Belczynski (2005) . As summarized in Yungelson (2005) , 
Tutukov & Yungelson quantified the fraction of all close 
intermediate-mass binaries that reach the various binary 
evolution endpoints, including SN-Ia. They found that the 
dominant channel for SNe-Ia is the DD merger channel, 
with 0.7% of all binaries with both ZAMS masses in the 
m = 0.8 — 11.8 range (the full range of masses that pro- 
duce WDs within a Hubble time) leading to a super-Chandra 
merger. Limiting the progenitor population to m = 3 — 8, 
the predicted fraction is fia = 0.14 (L. Yungelson, private 
communication), considerably less than the ~ 100% I have 
derived from the observations. The SD channel was found to 
be an order of magnitude less efflcient in producing SNe-Ia. 

Hachisu et al. (1996, 1999a,b; see also Kobayashi et 
al. 1998, 2000) have argued against DD mergers as SN-Ia 
progenitors, and have advocated, instead, two different, co- 
existing, SD paths to SN-Ia explosion: a primary C/O WD 
with ZAMS mass m = 5 — 8, accreting from a red-giant sec- 
ondary of initial mass m = 0.8 — 1.5;ora primary of ZAMS 
mass m = 5.5 — 8.5 that accretes from a ZAMS m = 1.8 — 3.4 
main sequence secondary0 Hachisu et al. claim that these 
SD channels are an order of magnitude more efficient than 
estimated by Yungelson & Livio (1998), i.e., comparable to 
the predictions for the DD channel favored by Tutukov & 
Yungelson (2002), given above. 

Han & Podsiadlowski (2004) have calculated more de- 
tailed numerical models for the WD + main sequence SD 
channel. They find a similar, though somewhat smaller, 
range of secondary initial masses, m ~ 2 — 3.4, that lead to 
an explosion. In terms of rates, they find that this SD chan- 
nel's efficiency is intermediate to that in the high results of 
Hachisu et al. (1999a) and the low results of Yungelson & 



^ Incidentally, Kobayashi et al. (2000) fit the element abundance 
patterns in the Solar neighborhood using the yields from their SD 
SN-Ia models and from core-collapse SNe. Their best fit value for 
the total ri, through both SD channels, is 7%, similar to the other 
estimates of rj, above. This can be considered a sixth independent 
estimate of rj, albeit a more model-dependent one. 



Livio (1998; see also Fedorova et al. 2004; see Han & Podsi- 
adlowski 2004, for a discussion of the sources of disagreement 
between the various studies.) The first SD channel, with a 
low-mass red giant secondary, was found to contribute neg- 
ligibly, in agreement with Yungelson & Livio (1998). 

In contrast. Hurley et al. (2002; see also Tout 2005) 
found very low SN-Ia rates from their DD channels but rates 
comparable to the higher rates cited above for some of their 
SD channels. These authors have also argued that acE, the 
parameter used to represent in a simplified way the effects 
of the physically complex common envelope stage of binary 
evolution, can assume values larger than unity. The con- 
stant QfCE is the fraction of the orbital energy that goes into 
driving off the envelope (see Nelemans & Tout 2005, for an 
alternative parametrized description of the common enve- 
lope stage). Most studies have assumed qce ~ 1, but these 
authors have calculated also models with qce ~ 3 which, 
they reason, are not unphysical, because additional sources 
of energy (thermal, from nuclear burning, or magnetic, from 
dynamo action) may be available, beyond the energy from 
rotation. Tout (2005) emphasizes that, by using high values 
of Qce, SN-Ia rates that are almost arbitrarily high can be 
obtained. 

To compare the observational estimates also to the SD 
theoretical predictions, we need to take into account the 
different progenitor mass ranges for this model - instead of 
m = 3 — 8 for both of the binary components, which I have 
used until now, we need to take for the primary m = 5 — 8, 
and for the secondary m — 0.8 — 3.4 ("generously" including 
also secondaries in the m — 1.5 — 2 range). Due to the smaller 
integration interval in the denominator, the rj estimates in 
Eqns. [T]|6] will grow by a factor 3.1, for the standard IMF 
slope. In other words, the fractions of all m — 5 — 8 stars that 
explode as SNe-Ia will be 3 times higher than the previous 
estimates. On the other hand, to describe the larger pool of 
secondaries, Eq. [7] needs to be replaced by 



J,^djV/dm)rfm/j^4,/-„ 
J^{dN/dm)dm 



f{q)dq 



(8) 



Compared to the previous result, fb ~ 1/3 — 1/6, this gives 
/b.SD ~ 1/2.5, weakly dependent on the power-law slope 
of /(g). Finally, because the mass ranges of the primaries 
and the secondaries are disjoint in this case, the duplicity 
factor, which we took as fd = 1/2 for the DD case, is now 1. 
Thus, among all the potential SD progenitor systems in the 
model of Hachisu et al. (1996, 1999a,b) - close binaries with 
a primary in the m = 5 — 8 range - the observed fraction 
of systems that produce SNe-Ia will still be 0.4-0.8 of the 
estimate obtained for the DD case, i.e., of order 100%, or 
even more. If we adopt the Han & Podsiadlowski (2004) 
conclusion, that only m = 2 — 3.4 secondaries can contribute 
non- negligibly to the SD channel, then /(,,sd ~ 1/4.3, and 
we come even closer to the DD result. 

Comparing the observational and theoretical results, 
the typical observed value of fia ~ 100% is an order mag- 
nitude larger than the DD prediction, and two orders larger 
than the the SD prediction of Tutukov & Yungelson (2002) . 
The SD rate predictions of Hachisu et al. ( 1999a, b) and 
Tout (2005, for acE = l) are comparable to the the DD 
rate prediction of Tutukov & Yungelson (2002), and the 
SD predictions of Han & Podsiadlowski (2004) are lower 
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than those of Hachisu et al. (1999a,b). Thus, all of these 
predictions are lower by an order of magnitude, or more, 
than the observed fraction. In fact, from fig. 6 of Han & 
Podsiadlowski (2004), one can deduc«[f| a predicted value of 
_B = (1 — 2) X 10"'* SNeM0~^, an order of magnitude lower 
than most of the measurements compiled in ^2.21 Even the 
lowest oberved value of 77 found above, ~ 2% (which is in- 
consistent with all other measurements), times the lowest 
values for {fafbfcfd)~^ ~ 12, still leads to fia ~ 24%, twice 
as large as the highest theoretical predictions. 

Previous comparisons of theory and observation, in the 
binary synthesis studies cited above, have focused on re- 
producing the Galactic SN-Ia rate, after assuming a Galac- 
tic star formation rate. Both of these rates are difficult 
to estimate, and are therefore highly uncertain. As a re- 
sult, a factor-few agreement was considered satisfactory by 
those studies. The modern measurements of extragalactic 
SN rates, and related observables, that I have used here 
make the order- of- magnitude discrepancy between theory 
and observations explicit. 

4.2 Possible solutions 

4-2.1 Observational errors 

To explain the discrepancy, we can re-examine some of the 
assumptions in either the observational derivation or in the 
theoretical prediction. First, we can consider the possibil- 
ity that the various observational estimates of SN-Ia rates 
have systematically erred on the high side. A SN-Ia rate 
may be systematically high due to contamination of a type- 
la sample by, e.g., core-collapse SNe or active galactic nu- 
clei. In fact, such contamination may be behind some of 
the discrepant SN-Ia rates measured at 2; ~ 0.7 by several 
groups (see discussion in Poznanski et al. 2007, and refer- 
ences therein), which lead to a fairly large range in the es- 
timates of riB- However, I have shown that even the lower 
estimates of rj are still high, compared with the theoretical 
expectations, and furthermore there is consistency with in- 
termediate values of rj ~ 15% from a variety of independent 
methods (not all of them involving SN rates). SN rates can 
also be overestimated if the sensitivity times of the surveys 
are underestimated, i.e., the surveys are more sensitive to 
faint SNe than they assume, which is quite unlikely to al- 
ways be the case. 

4-2.2 Mass-ratio distribution 

A flat distribution in the binary mass-ratio, f{q), was as- 
sumed for close binaries in the Tutukov & Yungelson (2002) 
simulations, while in the observational analysis I considered 
distributions that are either flat or mildly rising toward low 
q. An inverted mass ratio distribution that favors large (« 1) 
values of q would raise ft, and thus lower the "observed" //„ 
(in an extreme version, f{q) would be a 5 function at q — 1, 
i.e., all close binaries are "twins"). In the simulations, such 

Like all of the above binary synthesis works, Han &; Podsi- 
adlowski (2004) followed Iben & Tutukov (1984), in assuming a 
Galactic star-formation rate of 1 star yr~^ with m > 0.8. Assum- 
ing our standard IMF, this translates to 5.8M0yr~^, and gives 
the stated result for B. 



"twinness" would increase the rate of occurrence of super- 
Chandra mergers (although such equal-mass mergers would 
also be more prone to dynamical instability leading to core 
collapse, rather than a SN-Ia; Nomoto & Iben 1985). There 
have been some indications for a peak in f{q) at g = 1 for 
short-period binaries (e.g., Tokovinin 2000; Halbwachs et al. 
2003; Pinsonneault & Stanek 2006), but it has been counter- 
argued these are due to observational selection effects, and 
are not supported by other studies (see above) . Pinsonneault 
& Stanek (2006) find that 55% of a sample of detached bi- 
nary systems that they study is in a "twin" peak in the q 
distribution. However, their twins all have primaries in the 
m = 12 — 23 mass range, so even if this is not a selection 
effect, the relevance to the present study is not certain. 

4.2.3 Binary separation distribution 

A similar type of solution is to assume that there are more 
close binaries than indicated by modern studies, and this 
would again work toward making the observational and the- 
oretical estimates of fia more consistent with each other. 
Again, there is no solid observational evidence, or theoret- 
ical motivation, for such a deviation from the fairly well- 
established Oepik distribution. 

4.2.4 A "middle-heavy" IMF? 

An alternate approach to solve the problem may be to hy- 
pothesize a quite different form for the IMF. For exam- 
ple, a "middle-heavy" IMF with a "bump" at intermedi- 
ate masses could provide a larger pool of SN-Ia progenitors, 
and would lower the observational estimates of fia. Inter- 
estingly, Fardal et al. (2007) have recently proposed such 
a "paunchy" IMF, from a completely different motivation. 
They found that measurements of the extragalactic back- 
ground light (EBL) and of the cosmic star formation history, 
SFR{t), may overpredict the observed local stellar density, 
if one of the standard IMFs is assumed. As an example of 
a paunchy IMF that can solve this problem, they proposed 
an IMF with a slope a — —1.7 in the m = 0.5 to m = 4 
range, falling to a = —2.7 at higher masses. They noted 
the difficulties of measuring the IMF in this intermediate- 
mass region (Kroupa 2002), and the fact that Sirianni et 
al. (2000) had obtained, for the R136 cluster in the Large 
Magellanic Cloud (LMC), an IMF slope of q = -1.3 in 
the 1.3 < m < 2.1 range. The IMF proposed by Fardal et 
al. (2007), when normalized to match the Salpeter IMF at 
m = 1, has 1.5 times as many m = 3 — 8 stars as the Salpeter 
case, and rjA and r]Fe would be reduced by a factor 0.68. If 
normalized at m = 8, the reduction, which would apply to 
rjB, is only by a factor 0.96. The ratio of core-collapse to 
la progenitors is reduced by 0.74, compared to the Salpeter 
case, and this would be the effect on r;ia/cc a-nd on ?7abund- 
Thus, the effect of this particular paunchy IMF is overall 
moderate, and a more extreme variant is required to reduce 
rj more significantly. 

However, direct IMF measurements in star-forming re- 
gions generally do not show such a feature. For example, 
Preibisch et al. (2002) measured the IMF in the m = 0.1-20 
range for the Upper Scorpius OB association. They found 
the IMF to be broadly consistent, in the m > 1 range, with 
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the standard Scalo (1998) and Kroupa (2002) IMFs. Their 
best fit slope was a = —2.8 ± 0.5 in the m = 0.6 — 2 range, 
and a = —2.6 ± 0.3 in the m = 2 — 20 range. There have 
been other claims for top-heavy IMFs in the cores of dense 
young clusters in our Galaxy and in the LMC, only to be 
contradicted by other studies of the same objects. For ex- 
ample, Stolte et al. (2005) found for the Arches cluster an 
extremely top-heavy IMF, with a = —1.7 for m = 6 — 16, 
flattening to a ~ —1 below m = 6. Kim et al. (2006), in con- 
trast, deduced for the same chistcr an initial a of —2.0 to 
—2.1 in the m = 1.3 — 50 range, while Espinoza et al. (2007) 
find a = —2.4 ±0.1 in the m = 4 — 90 range, fully consistent 
with a Salpeter slope. Portogics Zwart et al. (2007) have re- 
cently used N-body simulations to argue that all of the data 
for the Arches cluster are consistent with a Salpeter IMF 
in the m = 1 — 100 range, but that the results for this and 
other young clusters are distorted by a combination of dy- 
namical and observational effects. Similar consistency with a 
Salpeter slope has been found in other clusters (e.g., Ascenso 
et al. 2007). In the Sag A* cluster in the Galactic center, 
Nayakshin & Sunyaev (2005) have used the low level of X- 
ray emission to deduce a factor-10 deficiency in the number 
of m < 3 stars, compared to the Orion Nebula cluster, and 
to expectations from a standard IMF. Near-infrared obser- 
vations by Maness et al. (2007) lead to similar conclusions. 
However, this case of top-heavy star formation is thought 
to be a peculiarity related to the proximity of the Sag A* 
cluster to the Galaxy's supcrmassivc black hole. Thus, it ap- 
pears that direct and persistent evidence for a non-standard 
IMF is scarce, and this is unlikely to be the solution the 
SN-Ia formation efficiency problem. 

4-2.5 Problems with the theoretical models 

Turning to the theoretical calculations, as noted above. Tout 
(2005) has made the point that values of acE sufficiently 
larger than unity can improve dramatically the efficiency of 
both the SD and DD channels. In his simulations, raising 
acE from 1 to 3 increases the SN-Ia rate by an order of 
magnitude, in some cases. If such an assumption is valid, 
this could be a solution to the problem. 

Finally, we may consider the more radical possibility 
that some of the more basic assumptions about SN-Ia pro- 
genitors are wrong. Both the DD and SD scenarios are not 
free of theoretical and observational problems. It has long 
been argued that the final stage of a DD merger will lead 
to an accretion-induced core collapse, rather than a SN-Ia 
(e.g., Nomoto & Iben 1985; Saio & Nomoto 2004, Guerrero 
et al. 2004), although there are opposing views that invoke 
rotation of the stellar surface to prevent this outcome (e.g., 
Piersanti et al. 2003). Au observational search for DD pro- 
genitor systems (Napiwotzki et al. 2004; Nelemans et al. 
2005) has turned up, among ~ 1000 WDs surveyed, few or 
no potential DD systems with a total mass exceeding Mch 
that will merge within a Hubble time. It is not yet clear 
if this low detection rate is significantly smaller than ex- 
pected if the DD scenario is to explain the Galactic SN-Ia 
rate (which is, itself, quite uncertain). For example, if the 
Galactic rate is lO^'^yr^^, there should be 10^ progenitor 
systems that will merge within 10^*^ yr. If there are 10^^ stars 
in the Galaxy, then, assuming a Kroupa (2001) IMF, of or- 
der 10%, or 10^°, will be WDs. Only one in a thousand WDs 



surveyed would therefore be expected to be a SN-la progen- 
itor, which is consistent, at this order-of-magnitude level, 
with the current survey. A more detailed calculation, tak- 
ing into account Galactic structure and star formation his- 
tory, survey completeness, and selection effects, is required 
in order to determine whether or not the observations are in 
conflict with the DD scenario. 

The SD scenario, in turn, has been criticized for its as- 
sumptions about the existence of ad hoc mechanisms that 
regulate the accretion flow on to the primary (e.g., Cas- 
sisi et al. 1998; Piersanti et al. 1999, 2000). Observation- 
ally, Badenes et al. (2007) have noted the absence, in seven 
nearby SN-Ia remnants, of the signatures of the strong wind 
from the accrctor that supposedly stabilises the accretion 
flow, and permits reaching Mch in the Hachisu et al. (1996) 
model. Prieto et al. (2007) have pointed out that a sub- 
tantial number of SNe-Ia have been discovered, by now, in 
low-mctallicity galaxies, in contrast to to the prediction by 
Kobayashi et al. (1998) that SNe-Ia cannot explode in such 
environments, due to a minimum metallicity that is required 
for the wind regulation mechanism to be effective. (Yoon et 
al. 2004 have proposed WD rotation as an alternative or 
additional mechanism for stabilising high accretion rates.) 
While evidence for circumstellar material, consistent with 
expectations from a wind from a red-giant companion, has 
been found in one recent normal SN-Ia (Patat et al. 2007), 
such material is not observed in another event (Simon et al. 
2007). For the remnant of Tycho's SN, which was a type- 
la (Badenes et al. 2006), there have been conflicting claims 
about the identification and nature of a remaining compan- 
ion star (Ruiz-Lapuente et al. 2004; Puhrmann 2005; Ihara 
et al. 2007). If, in the end, no companion is found, this would 
be another problem for the SD picture. 

4-2.6 Single-star SN-Ia progenitors? 

A speculative direction that may be worth considering is 
that, perhaps, SNe-Ia are not descended from binary sys- 
tems. Iben & Renzini (1983) have reviewed several stel- 
lar evolution studies of intermediate-mass stars, in which 
mass loss has not reduced the remaining stellar mass be- 
low ^ Mch at the time of carbon ignition in the core. Such 
single stars may undergo a thermonuclear runaway that dis- 
rupts the star, but Iben & Renzini (1983) concluded that, 
due to the many assumptions and approximations, "one is 
not left with an overwhelming sense of confidence in the 
detailed results". Such an explosion would be expected to 
take place still within a hydrogen-rich envelope, and hence 
would not resemble a normal SN-la. However, Tout (2005) 
has sketched the possibility that some not-fuUy understood 
aspect of stellar evolution leads all stars in some initial mass 
range to avoid mass loss on the red giant branch, but to then 
completely lose their hydrogen envelopes on the asymptotic 
giant branch (AGB), proceeding then to a Chandra-mass 
thermonuclear runaway of their cores in such a "single-star" 
SN-la channel. In a variant of this idea, Waldman, Yungel- 
son, & Barkat (2007) have postulated that interaction with 
a binary companion is responsible for the stripping of the 
envelope, but the stripped star then goes on to explode as 
a single star. 

Observationally, recent metal abundance measurements 
in Kepler's SN remnant strongly suggest a type-la explosion 
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(e.g., Blair et al. 2007; Reynolds et al. 2007). However, the 
~ 300 km velocity of the remnant as a whole away from 
the Galactic plane, together with evidence for a circumstel- 
lar medium from a massive progenitor, are difficult to under- 
stand in the context of a binary progenitor system (Reynolds 
et al. 2007). Aubourg et al. (2007) have modeled the stellar 
populations of SN-Ia host galaxies, and concluded that the 
prompt SN-la component explodes within 70 Myr of star for- 
mation, implying a primary initial mass range of m = 6 — 80 
For such a narrow and high range, all of the estimates of 
through the prompt channel will become even higher than 
found earlier. For example, rfia/cc will be 67 ± 15%. Even 
with the most conservative binary parameters (including a 
substantial "twin" population, which could raise /;,, but at 
the same time would restore the duplicity factor, fd = 1/2), 
it is then difficult to avoid an excess of SNe-Ia compared 
to progenitor binary systems. A single-star channel would 
solve this problem. 

A good place to test for the existence of such single- 
star progenitors could be the LMC. Taking a value of _B ~ 
1 X 10~^Mq ^ for the SN-Ia rate per star formation rate (see 
§2.2), which should be appropriate for star-forming dwarf 
galaxies like the LMC, and a SFR value of O.SMoyr"^ (Ken- 
nicutt et al. 1995), an LMC SN-Ia rate of ~ 3 x 10"*yr"^ is 
obtained. Another estimate can be made if the 2:1 to 4:1 ra- 
tio (Mannucci et al. 2005) of core-collapse to SN-Ia rates in 
star-forming galaxies holds for the LMC, then amongst the 
45 to 76 SN remnants in the LMC (Sasaki, Haberl, & Pietsch 
2000; Bojicic et al. 2007), roughly 9-25 should be the rem- 
nants of SNe-Ia. Divided by a typical 10* yr lifetime of a 
remnant, this gives a SN-Ia rate of ~ (9 — 25) x 10~*yr~^. 
The AGE stage of ZAMS m = 6 - 7 stars lasts of order 
10^ yr (e.g., Girardi et al. 2000). Thus, if SN-Ia progeni- 
tors are AGE stars that have lost their hydrogen envelopes, 
there could be a few tens or hundreds of such stars in the 
LMC. Calculated colors and spectra of such objects could 
be compared to observations by searching for candidate pro- 
genitors in photometric catalogues of luminous stars in the 
LMC (Cioni et al. 2000; Nikolaev & Weinberg 2000). 

To summarize, I have collected measurements perti- 
nent to five independent methods for estimating the frac- 
tion of potential SN-Ia progenitors that explode as SNe- 
Ia. I have compared them self-consistently, and studied the 
dependence of the fraction on the assumed IMF and on 
other parameters. Modern measurements of SN rates and 

* The SN sample of Aubourg et al. (2007) is a literature compi- 
lation from several SN surveys and from lAU circulars, for which 
selection effects and "control times" (the time during which a SN- 
Ia could have been detected in a given galaxy) are not known. In 
their derivation of a SN-Ia rate as a function of star formation 
rate, these authors implicitly assume that the control times and 
the foUowup efficiencies are the same for SNe in both the star- 
forming and the quiescent galaxies in the surveys that produced 
the SN sample. At some level, this assumption likely fails, and this 
casts some doubt on the accuracy of their result. For example, 
star- forming galaxies tend to host more luminous, more slowly 
fading, SNe-Ia, compared to early-type galaxies (see, e.g., Gal- 
lagher et al. 2005); such SNe are more likely to be detected, and 
more likely to be confirmed spectroscopically as type-la's, artifi- 
cially biasing high the derived SN-Ia rate in star-forming galaxies. 
If some of the surveys specifically targeted late-type galaxies, the 
bias would be even greater. 



related observables leave little doubt that a large fraction of 
the intermediate-mass stellar population explodes as SNe- 
Ia, larger than would be expected based on most SD and 
DD models. A critical reappraisal of our ideas about SN-Ia 
progenitors may be in order. 
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